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В данной работе был проведен анализ сигнала электрокардиограммы при помощи методов вейвлет- 
преобразования с целью выявления диагностических признаков и локальных особенностей сигнала. 


Введение 


На сегодняшний день одним из самых распространенных методов диагностики 
сердечно-сосудистых заболеваний является электрокардиография. 

Электрокардиограмма (ЭКГ) характеризуется набором зубцов, по временным и 
амплитудным параметрам которых ставится диагноз. Увеличение амплитуды зубцов 
ЭКГ говорит о гипертрофии соответствующих отделов сердца, которая наблюдается 
при некоторых заболеваниях сердца и при гипертонической болезни [1]. По ЭКГ также 
можно оценить источник (так называемый водитель) ритма, регулярность сердечных 
сокращений, их частоту — это имеет большое значение для диагностики различных 
аритмий. По продолжительности различных интервалов и зубцов ЭКГ можно судить 
об изменениях сердечной проводимости. Изменения конечной части желудочкового 
комплекса (интервал $Т и зубец Т) позволяют врачу определить наличие или отсутствие 
ишемических изменений в сердце (нарушение кровоснабжения). 

Выявление и оценка ключевых диагностических признаков электрокардиосигнала 
обеспечивает обнаружение нестабильности миокарда и диагностики угрожающих жизни 
нарушений ритма [2]. 


1 Постановка задачи 


Целью работы является создание программного обеспечения, позволяющего 
облегчить работу врача-функционалиста при постановке диагноза, а именно — при 
обработке электрокардиосигнала обнаружить и наглядно представить его ключевые 
диагностические признаки и выявить локальные особенности для определения наиболее 
вероятной патологии в работе сердца. 

Сигналы ЭКГ нестационарны и имеют сложные частотно-временные характе- 
ристики. Как правило, такие сигналы состоят из близких по времени, короткоживущих 
высокочастотных компонент и долговременных, близких по частоте низкочастотных 
компонент. 

Для анализа таких сигналов необходим метод, способный обеспечить хорошее 
разрешение и по частоте (для локализации низкочастотных составляющих), и по 
времени (для разрешения компонент высокой частоты). 
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Существует два подхода к анализу нестационарных сигналов такого типа. 

Первый — локальное преобразование Фурье (зво-Ите Еоицмег апп). 
Используя данный подход, мы работаем с нестационарным сигналом, как со стацио- 
нарным, предварительно разбив его на сегменты («окна»), статистика которых не 
меняется со временем. 

Второй подход — вейвлет-преобразование, которое представляет собой разложение 
нестационарного сигнала по набору базисных функций, определенных на интервале 
более коротком, чем продолжительность кардиосигнала. При этом все функции 
набора порождаются с помощью двухпараметрического преобразования (сдвига по оси 
времени и изменения масштаба) одной исходной функции (материнского вейвлета). 

Операция умножения на «окно» есть в самой базисной функции, при этом про- 
исходит адаптация «окна» к сигналу при изменении масштаба. 

Наиболее перспективным представляется использование вейвлет-преобразования 
при анализе острых состояний, когда ведется поиск резких изменений структуры сер- 
дечного ритма, а также для обнаружения непериодичных локальных низкоамплитудных 
особенностей сигнала ЭКГ [3]. 


2 Структура системы анализа ЭКГ 
на основе вейвлет-технологий 


Система (рис. 1) состоит из 4 основных блоков, каждый из которых реализован 
в соответствующей функции: 


Программа анализа ЭКГ в базисе вейвлетов 


Графический интерфейс пользователя 


Очистка сигнала от 
Загрузка сигнала искажений и 
неинформативных 
составляющих 


Реализация Реализация 


предварительного спектрального 
анализа ЭКГ анализа ЭКГ 


Рисунок 1 — Структура разработанной системы 


— блок «Загрузка сигнала» осуществляет загрузку та{-файла с данными сиг- 
нала ЭКГ; 

— блок «Очистка сигнала от искажений и неинформативных составляющих» 
позволяет из исходного сигнала удалить искажения, вызванные случайными помеха- 
ми при регистрации ЭКГ; 
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— блок «Реализация предварительного анализа ЭКГ» выполняет кратномас- 
штабный анализ (КМА) очищенного на предварительном этапе сигнала, а затем 
пользователю предоставляется возможность провести сравнение полученных резуль- 
татов с набором исходных показателей. Набор исходных показателей был получен в 
ходе исследования данных вейвлет-разложения для группы здоровых лиц и для 
группы лиц с инфарктом миокарда. Результатом работы этого блока является графи- 
ческое представление результатов КМА кардиосигнала данного пациента и набора 
исходных показателей на одной плоскости; 

— блок «Реализация спектрального анализа ЭКГ» выполняет процедуру непре- 
рывного вейвлет-преобразования сигнала при помощи вейвлетов Добеши-4 и Морле 
на заданных пользователем масштабах. Результатом работы этого блока является 
графическое представление результатов в виде спектрограммы. 


3 Экспериментальные результаты 


Был исследован сигнал ЭКГ в Ш отведении с частотой дискретизации 1000 Гц. 
Базовыми вейвлетами при обработке ЭКГ-сигналов были выбраны вейвлет Добеши-4 
(центральная частота — 0,7143 Гц) и вейвлет Морле (центральная частота — 0,8125Гиц). 

В работе представлены результаты обработки сигнала ЭКГ пациента с инфарктом 
миокарда. Исходный загруженный и очищенный от искажений сигнал имеет ярко 
выраженные отрицательные К-пики, что свойственно инфаркту миокарда. 

На этапе предварительного анализа ЭКГ выполняется кратномасштабный анализ 
заданного сигнала и производится отображение полученных данных на одной плоскости 
с параметрами двух исследованных ранее групп пациентов для дальнейшего визу- 
ального анализа (рис. 2). 
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Рисунок 2 — Результат выполнения предварительного анализа ЭКГ 


На рис. 2 достаточно точно отображается схожесть показателей вариабельности 
ЭКГ с показателями исследованных ранее групп пациентов, в частности (не зная 
диагноза) можно предположить, что у данного пациента имеются серьезные нарушения 
сердечной деятельности. Результаты предварительного анализа не являются материалом 
для диагностики, а являются как бы «подсказкой» и указывают врачу на общее 
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состояние сердечно-сосудистой системы человека в целом, а также указывают, на 
какие составляющие спектра необходимо обратить внимание. 

При реализации спектрального анализа ЭКГ выполняется процедура непрерывного 
вейвлет-преобразования сигнала при помощи вейвлета Добеши-4 и вейвлета Морле 
на заданных пользователем диапазонах частот. Результатом выполнения является изо- 
бражение (рис. 3), иллюстрирующее частотно-временные характеристики сигнала. 
По оси абсцисс откладывается время, по оси ординат — частота, а абсолютное значение 
вейвлет-преобразования для конкретной пары время-частота определяет цвет, которым 
данный результат будет отображен (чем в большей степени та или иная частота присут- 
ствует в сигнале в конкретный момент времени, тем светлее будет оттенок) [4], [5]. 

На спектрограмме можно выявить локальные особенности сигнала, соответству- 
ющие определенным патологиям сердечной деятельности. На малых масштабах вей- 
влетов (при широком диапазоне частот от 1 до 300 Гц) достаточно хорошо проявляется 
комплекс ОК$ (последовательность зубцов на кардиограмме, отображающая фазу со- 
кращения желудочков) и его динамика во времени. 
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Рисунок 3 — Результат выполнения спектрального анализа ЭКГ 


В диапазоне менее 3 Гц (при больших масштабах вейвлета) непрерывное вейвлет- 
преобразование (НВП) на основе вейвлета Морле не отражает особенности спектра 
сигнала, в отличие от НВП на основе вейвлета Добеши-4. Следовательно, для обработки 
ЭКГ-сигналов предпочтительнее использовать вейвлет Добеши-4, что позволит, не 
прибегая к пульсометрии, либо к дополнительной обработке ЭКГ, исследовать низко- 
амплитудные составляющие сигнала (рис. 4). 
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Рисунок 4 — Результат спектрального анализа ЭКГ при больших масштабах вейвлетов 
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Разработанный программный продукт обладает возможностью наращивания 
функциональности. Для более полного, всестороннего анализа ЭКГ и точной поста- 
новки диагноза данный программный продукт возможно дополнить блоками филь- 
трации и сжатия исходного сигнала, статистического, геометрического анализа 
кардиограмм, методом корреляционной ритмографии. Также важным направлением 
развития данного программного обеспечения является выявление и определение ко- 
личественных параметров, подтверждающих заключение врача, сделанное на основе 
визуальной интерпретации. 

На данном этапе разработанный программный продукт может быть использован 
в качестве дополнения и расширения уже имеющихся программно-аппаратных 
комплексов для анализа ЭКГ и диагностики сердечно-сосудистых заболеваний. 


Выводы 


Результатом проделанной работы является написанный на языке программирования 
пакета МАТГАВ 7.3.0 программный продукт, который реализует предварительный 
анализ ЭКГ и ее спектральный анализ. Полученные результаты позволяют обнаружить 
непериодичные локальные низкоамплитудные особенности сигнала ЭКГ, а также 
предоставляют возможности для исследования кардиограммы в диапазонах свыше 
100 Гц и менее 3 Гц с целью поиска новых диагностических признаков. 
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